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The experimental El strength distribution below 4 MeV in rare-earth nuclei suggests a local 
breaking of isospin symmetry. In addition to the octupole states, additional 1~ states with 

enhanced El strength have been observed in rare-earth nuclei by means of ( 7 , 7 ’) experiments. 
By reproducing the experimental results, the spdf interacting boson model calculations provide 
further evidence for the formation of an a cluster in medium-mass nuclei and might provide a new 
understanding of the origin of low-lying El strength. 

PACS numbers: 21.10.Re, 21.60.Ev, 21.60.Gx, 25.20.Dc 


The atomic nucleus is a unique laboratory which allows 
us to study how matter is built from the smallest scales 
(10“^®m) to stellar scales. In this mesoscopic system 
the symmetry between protons and neutrons, i.e. global 
isospin symmetry [i| , is a fundamental assumption. How¬ 
ever, already in 1985 one of us (F. lachello) proposed 
that even at lower excitation energies local rather than 
global isospin symmetry is realized and that enhanced El 
transitions could possibly further test the mechanisms 
by which this local symmetry is broken in atomic nu¬ 
clei Q. Below the particle-emission threshold, two kinds 
of El excitations have been intensive^ studied during 
the last decade, the octupole modes [1,I3 and the pygmy 
dipole resonance (PDR) [^, see Figs. □b) and [He). For 
both modes, it has been pointed out that they could have 
an important impact on fundamental physical properties, 
e.g., on enhancing the sensitivity for permanent electric 
dipole moment measurements in the case of the octupole 
mode J6| and on the equation of state in the case of the 
PDR 0. 

In both cases, the nonuniform distribution of protons 
and neutrons in simplified macroscopic models gives rise 
to enhanced El transitions at energies below the neu¬ 
tron separation energy. A clustering mode has also been 
proposed by one of us 0 where the a cluster might be 
considered as the simplest and energetically most favored 
realization, see Fig. □a). This mode would give rise to 
enhanced El transitions due to the oscillation of the a 
cluster relative to the remaining bulk 0,0 . 


In physics, clustering phenomena are of interest in 
many fields. For exarnple on the microscopic scale ul¬ 
tracold gases in traps Q, the electron-hole-exciton sys¬ 
tem in excited semiconductors fl^ . El or cluster systems 
for quantum computation 12|, [iSj are studied and on the 
macroscomc scale, e.g., clusters of stars [1^ and galax- 


16l | are of recent interest. In nuclear physics 


les 

a cluste ring is a well-established feature o f lig hter nu- 
[1^ [ij e.g., [l3-E| and EKIJ, and its 


clei 


implications for the El strength have been discussed [H- 


[20. First strong indications of a ^°®Pb-|-a system have 
been observed by means of enhanced El transitions be¬ 
tween excited states of ^^^Po [27| and also the possibility 
of a-cluster states in the actinides was discussed long 
ago, e.g.. Refs. 0 28-30. Receiitly, an exploratory cal¬ 
culation for ^^^Po was presented [3l|, indicating the exis¬ 
tence of ^°®Pb-|-a configurations when four-particle cor¬ 
relations are added to the shell-model calculations. This 
calculation provided a first hint at how to extend the well- 
established Tohsaki-Horiuchi-Schuck-Rdpke wave func¬ 
tion concept used for a-like condensates in light nu¬ 
clei 2lj, 13214341 to heavier nuclei. However, the general 
existence of a clustering in heavier nuclei remains an 
open question. Since a clustering in nuclei provides in¬ 
teresting insights about the formation of bosonic clusters 
in strongly coupled fermionic systems, indentifying new 
signatures of a clustering in heavier nuclei is of general 
scientific interest. 

In ^^'^Nd strong reduced a widths were observed for 
a group of states with J^= O’*", 2+, and 4+ between 2 
and 2.5 MeV in the {d,^Li) reaction, which was inter¬ 
preted as new evidence for a rotational band built upon 
an a-cluster state [30 . These results triggered systematic 
(7,7’) studies to find candidates for the expected J'"= 1“ 
states of the a-clustering mode 0. Indeed, several en¬ 
hanced El transitions were observed in 36- 

110, 148-154gj^ 156-160Q(J 

other rare-earth nuclei 40 below 4 MeV. First at- 



Figure 1. (Color online) Macroscopic interpretations of the 
different low-lying dipole modes. (a) a-clustering mode, 
(b) octupole mode, and (c) neutron-skin oscillation (pygmy 
dipole resonance (PDR)). 
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tempts to qualitatively describe the experimentally ob¬ 
served peculiar increase of the El strength connected 
with the proposed quadrupole-octupole coupled (QOC) 
Ij" state 5 towards shell closure s iS in a microscopic 
manner were presented in Refs. |44l - l50l| . There, it was 
proposed that one-particle-one-hole {Ip—lh) admixtures, 
originating from the tail of the electric giant dipole reso¬ 
nance (GDR), into the QOC state could explain the en¬ 
hanced El strength and that destructive interference of 
both components causes the decrease when moving away 
from the shell closure. Modern energy-density functional 
(EDF) calculations of El strength in transitional rare- 
earth nuclei were performed in Ref. (Hlj. However, the 
authors pointed out that the study of multiphonon states, 
e.g., 2+0 3“, with the EDF method is challenging and 
these have consequently only been studied in the spheri¬ 
cal semimagic N = 82 isotones using the second random- 
phase approximation (SRPA) [^, where they observed 
a complex microscopic structure of the 1 ]" state. 

It has been realized that cluster states form an addi¬ 
tional phase of nucleonic matter besides the mean-field- 
type states, see, e.g., Ref. and references therein. 
Recent theoretical studies in the framework of the EDF 
method have once more shown that different cluster con¬ 
figurations can appear as excited configurations in nuclei 
ranging from ®Be up to "^^Ca [s^. The algebraic clus¬ 
ter model [l^, [13 has been very successful to describe 
these configurations in light nuclei. Here, two-body clus¬ 
ters are described in terms of the algebra of 17(4). In 
fact, this corresponds to the sp version of the algebraic 
interacting boson model (IBM) and it has been 

proposed a long time ago that the p-boson is related to 
a-cluster configurations as shown in Fig. [I]a) [1, 281 and 
is expected to occur at the surface of heavier nuclei where 
the density is low, compare also Ref. 13ll. 

In the present study, the spdf IBM [5J-l5^ has been 
adopted to systematically study the low-lying J^= 1 “ 
states in the Nd isotopes and other rare-earth nuclei. 
This model has already been successfully applied to de¬ 
scribe El excitations related to octupole degrees of free¬ 
dom in the rare earths 55, 5^ and actinides The 

proximity of the additional and still not interpreted 1 “ 
states up to 4MeV to the octupole 1“ states indicates 
the need for full spdf IBM calculations. 

The following Hamiltonian has been chosen, which is 
a natural extension of the Hgd Hamiltonian [55| : 


Espdf — T Cpfip -t- efflf spdf ' Qspdf 

r / -N (3) / V (3)1 (0) 


03 




(1) 


with Cd, Cp, and e/ being the boson energies and fid, 
fip, and fif the number operators, respectively. The 
quadrupole interaction strength k of the quadrupole op¬ 
erator Qspdf is the same for describing positive- and 


negative-parity states simultaneously. To account for 
experimental anharmonicities, the 0(5) Casimir opera¬ 
tor [Slell has also been added. 03 corresponds to the 
strength of this 1= 3 interaction. We note explicitly that 
the choice of Hspdf describes sd and pf states separately 
and does not mix positive- and negative-parity boson 
states; i.e. besides the negative-parity boson energies all 
parameters are determined according to the signatures 
of well-established positive-parity collective states. For 
positive-parity states, these included the energies of the 
2f, 4+, 0^, and 2+ states as well as their reduced tran¬ 
sition strengths and 7 -decay branching ratios. For this 
approach see, e.g.. Ref. [^. The p- and /-boson ener¬ 
gies were fixed to describe the energies of the 1/ and 3/ 
states. To describe El transitions, the full one-body El 
operator was adopted: 


f{El) =ei[xsp(s'*'/-bp'''s)^^^ -b (p'*'cJ- 
+ Xdf{d'^f + Pd)’'^^], 




( 2 ) 


and its parameters were smoothly varied, see Fig. m{c), 
to reproduce the reduced transition strengths related to 
the decays of the 1/ and 3/ states to the yrast positive- 
parity states. 

The calculations for the Nd isotopes are compared to 
the experimental results in Fig. [2] and Table HI Both the 
evolution of the excitation energy of the 1 / state as well 
as the B{E1) / strength are nicely described. Typically, 
this state has been interpreted as the K= 0 projection 
of the one-octupole phonon excitation in deformed nu¬ 
clei and as a candidate for the J^= 1 “ member of the 
{2^ 0 3/) quintuplet in vibrational and spherical nu¬ 
clei Q. Indeed, when identifying the basis states by 
means of |[ns][np][nd][n/]) [^, then the 1 / in ^“^^Nd is 
found to be dominated with 55 % by the |40I1) configu¬ 
ration. This corresponds to the two-phonon 1“ state [65j. 
In addition, several J= 1 states are observed up to 
~ 4 MeV in experiment, for which negative parity has 
been assigned either by means of their K quantum num¬ 
ber assignment based on their 7 -decay properties Idol, i.e . 
AK— 0, or by direct parity measurements [s^, 1^1^. 
These states are marked black in the top panel of Fig. [2] 
For the other J= 1 states, the experimental situation 
without parity assignment is more complex since also the 
scissors mode contributes to the dipole distribution in de¬ 
formed nuclei at E^ ~ 3 MeV . Based on the present 
experimental data, it cannot be excluded that some of 
the J= 1 states do have positive parity. These states 
are marked in red. Furthermore, J^= 1~ candidates ob¬ 
served in inelastic scattering experiments are marked in 
gray [ 6 ^. Their scattering cross section has been scaled 
with respect to the observed one for the 1 / state and 
its respective B{E1) value. Most likely, they have not 
been observed in the ( 7 , 7 ') experiments due to the high 
background at low energies usually present in NRF exper- 
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Figure 2. (Color online) El distribution in To p p anel: Experimentally firmly assigned = 1“ states or assigned 

AK= 0 ground-state transitions are marked in black |3a - [3a . States with J=1 assignment but no parity and no AK=0 
assignment are marked in red [3^ . Dipole states measured in Ref. by means of {p,p') and (d, d') reactions are presented in 
gray. Midpanel: B{E1) strengths and bottom panel: rip/rif ratios predicted by the spdf IBM calculations, respectively. 


iments with bremsstrahlung. Keeping these limitations 
in mind, a good agreement between experiment and the 
present calculations is recognized for the additional El 
strength. 


Table I. Excitation energy and El strength for the 1)" state 
as well as centroid energy and summed El strength for the 
remaining 1“ states in the Nd isotopes. 


A 

[MeV] 

B(E1)^- t 

[10-® e^fm^] 

E,- 

u,/a Ij 

[MeV] 

[10-® e^m^j 

142 

3.4 

Experiment 

17(4) 


144 

2.19 

6 .1(10) 

3.5 

19.5(6)^ 

146 

1.38 

5(2) 

3.3 

14(2) 

148 

1.02 

14(6) 

3.3 

30(2)^ 

150 

0.85 

15(7) 

3.2 

18.3(11)^ 

142 

3.3 

spdf IBM 

19 


144 

2.00 

6.7 

3.4 

11.5 

146 

1.38 

5 

3.4 

12 

148 

1.02 

9 

3.3 

19 

150 

0.85 

19 

3.3 

11.2 


^ 10.4(5) for firm J'^ = 1~ states and AK = 0 assignment. 
^ 17.6(13) for J = 1 states with K assignment. 

^ 11.4(8) for firm = 1~ states and AK = 0 assignment. 


To further study the structure of the states, it is neces¬ 
sary to look beyond the strength distribution. The lowest 
panel of Fig. Oshows the p- and /-boson fraction npjnf. 
For illustrative purposes, only states are shown having a 
B{E1) larger than 0.1 x 10“^ e^fm^. In all Nd isotopes, 
states are observed which show an enhanced p-boson con¬ 
tent, i.e. npjnf > 1, and which make up for a great part 
of the enhanced El strength. 

Another remarkable observation is the stability of the 
centroid energy and summed strength of this additional 


El strength which is located around 3.4 MeV, see TableHl 
This automatically raises the question if and how this 
strength is connected to the proposed (2+(g)3“)i- state in 
i^^Nd at 3.4 MeV with B(F;1) t= 17(4) ■ e^fm^ [67 1. 
This structure assignment was based upon two facts: 
(1) E^^ « E 2 + E E 3 - and (2) B{E2;1- 3/) « 

B[E2\2jj Oj') [67|, supposed key signatures of a QOC 
state. However, a final assignment also requires another 
proof, namely, that B{E3-, 1“ —>■ 2)*") ss B{Ei;3j O)*"). 
The inspection of Fig.[2]shows that the 1/ state in ^"^^Nd 
has a non-neglibile p-boson content (np/n/=0.8). The 
|5100) configuration accounts for 28% of the total wave 
function (<10% in ^^^Nd), while |4011) contributes with 
36%. 

Furthermore, the (s^p -I- contribution to the 



Figure 3. (Color online) El matrix element (a) 
(ir||f’(Al)||0/)2 and (b) (3r||f'(Al)||2+)2 for 
Shown are the bare contributions to the spdf IBM El op¬ 
erator of Eq. [ 2 I No fit parameters have been adjusted, (c) 
Evolution of the fit parameters Xap and Xdf in Eq.[2]used to 
describe the experimental £11 strength. 
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Figure 4. (Color online) Same as Fig. [2]but for and In addition, states with an experimentally 

ambiguous 7 -decay branching are marked in yellow. 


(1)"||T(£'1)||0)'')^ matrix element in Fig. |3](a) clearly re¬ 
sembles the experimentally determined parabolic behav¬ 
ior of the El strength, see, e.g., Ref. [43| and references 
therein. Realizing the sp-boson space as the bosonic 
manifestation of Lf(4), which was proposed to describe 
two-body clusters [ISl, l22| , no complex mixing with GDR 
components is needed. As expected from macroscopic 
considerations, the other matrix elements nearly vanish 
at shell closure and cannot explain the El strength in¬ 
crease at shell closure. The p-boson structure also ex¬ 
plains the near equality of the two mentioned B{E2) val¬ 
ues since they correspond to E2 transitions of the type 
s^d and /^p, which lead to comparable strength when 
considering the spdf-IBM E2 operator . We also note 
that comparable El matrix elements, found in this Let¬ 
ter, can explain the suggested empirical correlation be¬ 
tween the B{E1; 1)" 0)^) and B{El;3i —>■ 2'^) values 

in vibrational nuclei [b^, see Fig.|3](a) and (b). 

To show the general importance of the p boson and 
the appearance of p-boson states all over the rare-earth 
region, we have also studied 

I62,i64j7y parameters for positive-parity states of 

the Dy isotopes were taken from Ref. [62|. The results 
are shown in Fig. SI As in the Nd isotopes, the general 
agreement with the experimentally observed El distri¬ 
bution is good. The similarities between and 

i46,i48]\Jj^ are in fact striking. Very strong p-boson states 
are observed in the Dy isotopes. 

The a-clustering mode, i.e. the p-boson mode, corre¬ 
sponds to 4-QP configurations, which were also proposed 
in the framework of the EDF a ppr oach to obtain a bet¬ 
ter description of the 1~ states Calculations of this 
kind have not been performed yet but the importance 
of the particle-particle interaction for semimagic nuclei 
has been pointed out in Ref. [i^. In Ref. [b^ it was 
shown that 2-QP components had to be admixed to the 
1^ states to describe the El strength, even though de¬ 
viations were found for the Nd and Sm isotopes. It was 
speculated that these admixtures could be mimicked by 


the p-boson. The clear evolution of from 1.5 MeV in 
^®°Nd to 4 MeV in ^'^^Nd and, especially, the independent 
evolution of the (1)"||T(£'1)||0)*')^ and (3)'||T(i?l)||2)^)^ 
matrix elements as well as the non-negligible p-boson 
contribution to the norm of the wave function, in con¬ 
trast to Ref. [b^ , are more likely to point to an evolution 
of a configuration complementary to the mean-field-like 
states and might provide a new understanding of the ori¬ 
gin of low-lying El strength in nuclei. 

In conclusion, the experimental El strength distribu¬ 
tion below 4 MeV suggests the presence of a new col¬ 
lective dipole mode different from the octupole mode in 
rare-earth nuclei ( 4 ^ , i.e. the a-clustering mode occuring 
at the surface of nuclei just above magic numbers ii. 
In this Letter, we have presented, for the first time, a 
systematic study of the El strength observed in rare- 
earth nuclei using the spdf IBM. In agreement with ex¬ 
periment, several enhanced El transitions are observed. 
The new class of excitations is closely related to the p 
boson. Furthermore, an alternative interpretation of the 
peculiar B{E1) increase at shell closure might be pos¬ 
sible in terms of the fundamental p-boson mode, i.e., a 
clustering [1, [ 2 ^ . Our studies suggest a close connection 
between El strength and a clustering Q, i.e., an addi¬ 
tional isospin-symmetry breaking component, in heavier 
nuclei and might hint at the general occurrence of this 
mode in nuclei. Microscopic calculations looking into the 
details of the wave functions as well as experiments sen¬ 
sitive to a structures, e.g., {d,^Li), (^Li,d), ( 7 , 7 ), and 
(e, e') reactions, are asked for. Furthermore, the ( 7 , 7 ') 
measurements should be extended up to the particle- 
emission threshold to study the complete low-lying El 
strength with evolving deformation. 

The experimental ( 7 , 7 ') data have been obtained in a 
very close collaboration of the Universities of Stuttgart, 
Darmstadt, and Koln. We thank all colleagues from this 
collaboration. Discussions with V. Derya, A. Hennig, and 
P. von Brentano are highly appreciated. This work was 
supported by the DFG (ZI 510/4-2) and in part by the 
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